Experimental Observation of Dramatic Differences in the 
Dynamic Response of Newtonian and Maxwellian Fluids 
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An experimental study of the dynamic response of a Newtonian fluid and a Maxwellian fluid 
under an oscillating pressure gradient is presented. Laser Doppler anemometry is used in order to 
determine the velocity of the fluid inside a cylindrical tube. In the case of the Newtonian fluid, the 
dissipative nature is observed. In the dynamic response of the Maxwellian fluid an enhancement at 
the frequencies predicted by theory is observed. 



PACS numbers; 47.60.+i, 47.50.+d, 83.60.Bc 



I. INTRODUCTION 



Since the concept of dynamic permeability was pro- 
posed several years ago P, 0, it has been extensively 
used in the study of practical problems such as petroleum 
recovery , soil-ground water transport ^ , fluids flow- 
ing in porous media 0| , acoustic waves in porous me- 
dia 7] , wave propagation in foams Q , fluid circulation in 
biological systems Q, etc. The use of the dynamic per- 
meability as originally proposed [3| however, is not quite 
correct for some of these systems since the fluid does 
not behave as a Newtonian fluid but shows viscoelas- 
tic features. Theoretical analyzes of viscoelastic fluids 
have recently predicted an interesting enhancement in the 
dynamic permeability of several orders of magnitude in 
comparison with the static permeability [l^, fill Il2 . [l3| . 
This behavior is due to the coupling between the elas- 
tic behavior of the fluid and the geometry of the con- 
tainer and is completely different from the pure dissipa- 
tive behavior of Newtonian fluids. An increase in the 
flow rate of a viscoelastic fluid flowing in a tube un- 
der oscillating conditions was discovered several years 
ago 0,0,0,01 but up to now, the experimental mea- 
surements in terms of a frequency-dependent response 
had not been performed. Some theoretical studies have 
explored interesting consequences of the enhancement of 
the dynamic permeability in different systems including 
viscoelastic fluids [llllllllliillSlillllllii. The 
variety of systems where an oscillating pressure gradi- 
ent can be imposed to a viscoelastic fluid is wide, for 
instance, from the oil recovery problem |24j to the dy- 
namic analysis of fluid transport in animals (including 
the human body H^). All previous theoretical deriva- 
tions of the dynamic permeability for viscoelastic fluids 
have been made in terms of the averaged velocity. In this 
work however, we perform a detailed analysis by measur- 
ing the dynamic response at the center of a fluid column 
that is moved in an oscillatory way and its local veloc- 
ity '2^ is measured using a laser Doppler anemometer 
(LDA). This allows us to test a local, simple model for 
both Newtonian and Maxwellian fluids and our results 
corroborate that the simple linear model can capture the 



differences between the dissipative (Newtonian fluid) and 
resonant (Maxwellian fluid) behaviors. 



II. LINEAR THEORY 



In the following, we present the pertinent equations 
used to obtain a local expression for the Maxwellian fluid 
response based on the results presented elsewhere |12j. 
We begin with the linearized momentum equation 
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and the continuity equation 

V-v = 0, 



(1) 



(2) 



for an incompressible fluid. In the above equations, p 
denotes the mass density of the fluid, v the fluid veloc- 
ity, p the pressure, and V* represents the viscous stress 
tensor. To ensure the validity of Eq. ^ it is necessary to 
consider fluids with low Reynolds numbers, in our exper- 
imental case: Re < 7x10"'^. The constitutive equation 
of the fluid that we use is the linearized Maxwell model: 
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where rj denotes the dynamic viscosity and t„j the re- 
laxation time. It is necessary to stress that the lin- 
earized Maxwell model constitutes a simplification re- 
quired to obtain simple analytic results that implies ne- 
glecting non-linear terms which may be important under 
certain circumstances. 

Manipulating the three previous equations, applying 
the Fourier transform, and using cylindrical coordinates, 
we obtain 
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der radius, P and V are the pressure and the velocity in 
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the frequency domain, respectively, lo is the angular fre- 
quency, and Jq is the cylindrical Bessel function of zeroth 
order. A more detailed derivation of these expressions 
can be found in 0, • 

Previous theoretical analyzes of the dynamic perme- 
ability of viscoelastic fluids deal with the above equations 
using an average velocity for the fluid. In what follows, 
we shall analyze the flow from a different point of view, 
one which is based on a "local response" defined as the 
response of the fluid to a pressure gradient in a particu- 
lar position inside the tube. This is required by the LDA 
that essentially measures quasi-local particle velocities. 
Thus, in order to analyze the local dynamic response, 
and to facilitate the measurements, we evaluate Eq. Q 
at the center of the cylinder (r = in cylindrical coordi- 
nates) where the amplitude of the velocity is maximum: 



(1 - ^CJ^„,) 



1 



Jo (/3a) 



dz 



The local dynamic response is thus defined as: 

V 
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where we have followed Darcy's generalized equation |lCll |. 
Substituting the velocity given by Eq. (|SJ), we obtain: 
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(7) 



From this expression one can recover the Newtonian be- 
havior with the substitution: tm — 0. In order to compare 
the last expression with experimental results, we need an 
expression for the pressure gradient. In our experimen- 
tal case, we have a harmonic oscillating column of fluid 
being moved by a piston coupled to a motor with ad- 
justable frequency. Thus, the pressure gradient can be 
easily described by: 



dpjt) 
dz 
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(8) 



which represents an oscillatory movement in the z direc- 
tion with a displacement amplitude equal to zq. We shall 
also measure the root mean square velocity, Vrmsi since 
the range of frequencies that can be resolved in this case 
is much wider than the range given by the Fourier trans- 
form of the velocity data. The experimental dynamic 
response of a fluid, (^™p(w)), is thus given by: 
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(9) 



Now, since the experimental velocity is harmonic, we 
can directly compare the value of the local dynamic re- 
sponse obtained using the previous expression for each 
experiment with the theoretical prediction from Eq. |(7J| . 
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FIG. 1: Schematic view of the setup used for the study of the 
dynamic response with an LDA measuring system 



III. EXPERIMENTAL SYSTEM 

The experimental device shown in Fig. 1, consists of a 
vertical cylinder filled with a fluid. The oscillating move- 
ment required to have a harmonic oscillating pressure 
gradient, c. /. Eq. JSJ, was produced with a piston that 
closes the base of the cylinder and is driven by a motor 
of variable frequency; as previously mentioned, the lo- 
cal velocity of the fluid is measured with a laser Doppler 
anemometer. 

The cylindrical tube, made of transparent acrylic, has 
an inner diameter of 5 cm and a length of 50 cm. In its 
lower end, this tube is joined by a clamp to a stainless 
steel piston skirt of the same inner diameter. Within the 
stainless steel skirt, there is a Teflon piston that is moved 
by a Siemens motor of 1 Hp, regulated by a Varispeed 
606PC3 system. This allows for a frequency of oscilla- 
tion that can be varied from 1.5 to 200 Hz. The acrylic 
cylinder over the skirt is contained in a square recipient 
with parallel walls of transparent acrylic and filled with 
glycerol having a refractive index similar to the one of 
the acrylic; this is in order to avoid the cylindrical aber- 
ration due to the fact that the cylinder acts as a lens. In 
this and other cases, when there are cylindrical configu- 
rations, these aberrations could cause a null intersection 
between the LDA beams or could result in the backscat- 
tering signal not being received by the photomultiplicr. 
The laser Doppler anemometer (LDA) used for the lo- 
cal velocity measurements is a widely known system . 
Our particular LDA consists of a probe (Dantec Fiber- 
Flow 60x20), a photomultiplicr, a Burst Spectrum Ana- 
lyzer {Dantec BSA 57N11), and an argon laser (Spectra 
Physics 177) emitting in four different wavelengths; only 
the 488 nm line (200 mW) is used, and the data are pro- 
cessed in a PC with the Dantec FlowManager software. 
The measurement volume is of 0.64 x 0.075 x 0.075 mm 
located at the center of the acrylic cylinder. The LDA 
probe was set into a positioning system (precision: 1 /im) 
placed on a different table from the rest of the setup 
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to prevent oscillations that could alter the experimental 
data. The rms-velocity is calculated from the collected 
data. 

Commercial glycerol was used as the Newtonian fluid 
to be tested; it has a viscosity of 1 Pas and a den- 
sity of 1250 kg/m^. The well-known cetylpyridinium 
chloride and sodium salicylate solution (CPyCl/NaSal, 
60/100) (2S . >29i] was used as the Maxwellian counter- 
part; its properties are p = 1050 kg/m^, rj = 60 Pas, 
and = 1.9 s |30l |. All the measurements were made 
within the 25±0.5°C interval in order to keep constant 
the properties of CPyCl/NaSal solution [s^]. The prop- 
erties of this solution and the diameter of the tube give a 
Deborah's number {a = pa^ / rjtm) of a = 0.0058. This 
value is much smaller than the critical value a < 11.64 
for the appearance of resonances predicted by the theo- 
retical model 0, |23| . The experimental frequency range 
is kept within the low Reynolds number regime by setting 
the piston movement amplitude to 0.8±0.05 mm (which 
assures that Re < 7 x 10^''). Under these circumstances 
we checked that the system can be described by a lin- 
earized balance momentum equation. Dantec 20 pm 
polyamid spheres were used as seeding particles; they 
remain suspended for a long time and cause a minimal 
standard deviation in the velocity data since the size of 
the particles is uniform. We can assure that the parti- 
cles follow the flow because, using Stokes's law (31j, the 
obtained limit frequency 



for which the particles follow an oscillation with a devi- 
ation up to 1 %, in amplitude, is much greater than the 
experimental frequencies (between 1.5 Hz and 15 Hz); 
/ = 9.7x10^ Hz for glycerol and / = 5.8x10^ Hz for the 
CPyCl/NaSal solution. The particle density, pp, is 10.3 
kg/m"^ and the radius of a particle, i?, is lxlO~^ m [3^ . 

To avoid transient perturbations, measurements were 
taken approximately 5 minutes after every change in 
the frequency. 

IV. RESULTS 

For each frequency 2, 000 velocity data were acquired 
with the LDA at the cylinder axis and 40 cm above 
the piston, where edge effects are negligible. Experi- 
ments were carried out in the oscillating frequency range 
[1.5, 15] ±0.1 Hz. In this frequency range the viscoelastic 
shear modulus is nearly real- valued (lossy elastic solid). 
The calculated rms-velocity from the LDA measurements 
for CPyCl/NaSal and for glycerol are shown in Fig 2, 
where one observes that the measured Maxwellian veloc- 
ity is dramatically different from the Newtonian velocity 
in the studied range. Clearly the rms-velocity increases 
with the frequency in both cases, but the velocity of the 
viscoelastic fluid shows well-defined peaks. An impor- 
tant point to be stressed is the fact that the velocity in 
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FIG. 2: Local ttos- velocity at the center of the cylinder. The 
circles correspond to the viscoelastic fluid and the triangles 
to the Newtonian fluid 



the glycerol case is a linear function with respect to the 
frequency which means a good quality of the mechani- 
cal piston movement. Moreover, this confirms that the 
movement of the piston is transmitted directly to the 
fluid. 

Originally inferred by the pioneers in the analysis of 
the dynamic permeability P, the dissipative nature 
of the glycerol response is clearly seen in Fig. 3, where 
the glycerol dynamic response and the frequency have 
been scaled, the first one by the value of the dynamic 
response at uj — Q and the second by the viscous time 
T — a} plr\. Also shown in Fig. 3 are the theoretical 
predictions derived from Eq. (7) (magnitude of ^ with 

tm = 0). 

Fig. 4 contains the dynamic response of the 
CPyCl/NaSal solution. Experimental values are plotted 
together with theoretical predictions (Eq.[7|). A dramatic 
change in the behavior of the dimensionless Maxwellian 
response with respect to the dissipative Newtonian case 
is observed, namely resonant frequencies appear. The 
agreement with Eq. is manifest: the positions of the 
peaks in Eq. |(7J) and in the experimental results occur 
at approximately the same values of the frequency cj, 
even though the amplitudes differ. We need to recall 
that the theory used for the prediction is a linear ap- 
proximation and neglects the convective term and the 
possible non linearities in the viscoelastic behavior of the 
fluid. Although we have selected the viscoelastic fluid 
in order to maintain the physical conditions where res- 
onance behavior appears, we were unable to assure the 
linear characteristic of the viscoelastic fluid. According 
to (33 the CPCl/NaSal solution behaves as a Maxwellian 
fluid if the shear rate is 7 < 0.6 s~^. Transforming Eq. 
4 to the time space and taking the radial derivative we 
have calculated the shear rate that depends on the radial 
coordinate and time. The maximum local shear rate in- 
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FIG. 3: Dynamic response for glycerol. Experimental values 
are shown by triangles (Eq. 9) and the line represents the 
theoretical prediction (magnitude of ^ from Eq. 7 with tm. = 
0) 



creases with frequency, e. g. ^2Hz < 1-5 s~^, jq.shz < 21 
s~^, 7ioffz < 70 s~^, where the subscripts indicate the 
frequency. Then, one of the possible causes of the dis- 
agreement in the response amplitude is because the fluid 
was locally under higher shear rates than its Maxwellian 
limit. This is in agreement with the fact that if we in- 
crease the frequency, then we find larger differences be- 
tween theory and experiment. Another reason for theo- 
retical and experimental differences could be compress- 
ible effects in the viscoelastic flow that also increase with 
frequency. 

These results constitute the first experimental evidence 
of the resonant behavior of a viscoelastic fluid and shows 
the validity of the theory p^Ell. 



ing the resonant behavior of the dynamic response of a 
Maxwellian fluid. The resonances appearing in the vis- 
coelastic response, as expected, are significantly higher 
than the monotonic decay in the Newtonian response. In 
both cases, the qualitatively agreement between theoret- 
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FIG. 4: Dynamic response for the CPyCl/NaSal solution. Ex- 
perimental values are shown by circles (Eq. 9) and the line 
represents the theoretical prediction (magnitude of ^ from Eq. 
7) 



ical predictions and the experimental results presented 
here is evident; the quantitative differences are probably 
due to non-linear effects of the viscoelastic fluid used. 
Our experimental results show that the dynamic response 
of viscoelastic fluids clearly exceeds the static response 
for some specific frequencies, in fair agreement with the 
theory for such systems. Finally, we want to stress that 
the results obtained in this work have a wide spectrum 
of applications, from oil recovery problems to the 
dynamical analysis of biological fluids ,3^ . 



V. CONCLUSIONS 

The local dynamic response of a fluid (either Newto- 
nian or Maxwellian) has been developed following pre- 
vious work P, 12, 0, ^3- This theoretical expression 
is compared with experimental LDA measurements in a 
harmonically-oscillating pressure system. The expected 
dissipative behavior of Newtonian fluids is confirmed. 
On the other hand and although the theoretical ap- 
proximation presented is a drastic simplification of the 
real viscoelastic behavior, it is still capable of reproduc- 
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